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Abstract. An explicit expression for the quadratic density- response function of a 
many-electron system is obtained in the framework of the time-dependent density- 
functional theory, in terms of the linear and quadratic density-response functions of 
noninteracting Kohn-Sham electrons and functional derivatives of the time-dependent 
exchange-correlation potential. This is used to evaluate the quadratic stopping power 
of a homogeneous electron gas for slow ions, which is demonstrated to be equivalent to 
that obtained up to second order in the ion charge in the framework of a fully nonlinear 
scattering approach. Numerical calculations are reported, thereby exploring the range 
of validity of quadratic-response theory. 
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1. Introduction 

The inelastic interaction of charged particles with matter is one of the fundamental 
problems of contemporary physics. It encompasses such phenomena as the stopping 
power of solids for moving ions, electron and positron energy-loss spectroscopy, inelastic 
low-energy electron diffraction, and hot-electron dynamics E] . 

A fruitful approach to the theoretical treatment of particle-solid interactions has 
proven to be the use of perturbation series expansions in powers of the projectile- 
target Coulomb interaction. For the description of many-electron targets, one typically 
introduces linear and quadratic density-response functions, which describe the electron 
density induced by external perturbations. 

While linear-response theory has proven successful in the description of the 
interaction of fast projectiles with solids, in the case of low projectile velocities and 
low electron densities a nonlinear description becomes quantitatively necessary [HIIUE]- 
Besides, there exist phenomena which cannot be explained in the framework of linear- 
response theory, an example being the existing difference between the scatterings of 
positively and negatively charged particles [El El- 

The random-phase approximation (RPA) has served as the natural starting point 
for the calculation of both linear [S] and quadratic |U E] density-response functions of 
the homogeneous electron gas. However, exchange and correlation (xc) effects, which 
are absent in the RPA, are known to be important for metallic electron densities j^j. 
The purpose of this paper is to derive in the framework of the time-dependent density- 
functional theory (TDDFT) [TTH ITT] an explicit expression for the quadratic density- 
response function of a many-electron system, which will then be used to evaluate the 
second-order energy loss per unit path length of charged particles moving through solid 
targets, i.e. the so-called stopping power of the target. 

Another approach to evaluate the energy loss of slow ions moving in a many-electron 
system is based on the ordinary formulation of scattering theory. In this approach 
^2jE3E], the stopping power for a heavy particle is determined in the low- velocity limit 
from the knowledge of the scattering phase-shifts, which can be obtained from a static 
nonlinearly screened potential by solving self-consistent ly the Kohn-Sham equation of 
density-functional theory (DFT) ^H]- Since these nonperturbative calculations include 
all orders in the projectile-target interaction, they represent an important standard to 
investigate the range of validity of perturbative expansions. Nonetheless, they have the 
limitation of being restricted to low velocities (v << vp, vp being the Fermi velocity) 
of recoilless probe particles moving in bulk materials. $ 

The interrelation of the perturbative-response and nonperturbative-scattering 
approaches in their overlapping range of applicability (low-velocity limit and small 
projectile charge) is both an interesting and non-trivial problem. The starting points of 

| An extension of DFT-based potential-scattering calculations to finite (although still small) projectile 
velocities has been reported in Ref. . The interrelation of this approach and the quadratic response 
theory at finite velocities is, however, beyond the scope of the present work. 
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these two schemes are completely different and there are no grounds to a priori assume 
equivalence between them. In this paper, we demonstrate that in the low-velocity limit 
and to second order in the external perturbation our quadratic-response formalism and 
the scattering approach are equivalent, thereby extending the RPA-based proof reported 
in Ref. [Sj to the general case where the xc effects are included. 

This paper is organized as follows. In Sec. 121 we derive in the framework of TDDFT 
a formally exact explicit expression for the quadratic density-response function of a 
many-electron system, in terms of the noninteracting Kohn-Sham linear and quadratic 
density-response functions and functional derivatives of the time-dependent xc potential. 
In Sec. El we derive basic expressions for the stopping power of a uniform electron gas, 
in the framework of both quadratic-response and nonperturbative-scattering schemes. 
The results of numerical calculations are presented in Sec. |H We use atomic units 
throughout, i.e., e 2 = h = m e = 1. 

2. Quadratic density response 

In the framework of TDDFT, Petersilka et al. J7j demonstrated that within linear- 
response theory the electron density n\ (r, t) induced in an arbitrary interacting many- 
electron system by the time-dependent external potential <ft ext (r, t) coincides with the 
electron density induced in the corresponding system of noninteracting Kohn-Sham 
electrons by the time-dependent effective potential 



where v(r, r') = l/|r — r'| is the bare Coulomb potential, and f xc [no\(r,t;r' ,t') is the 
functional derivative of the time-dependent xc potential \4 c [n](r,t) of TDDFT, to be 
evaluated at the unperturbed static electron density no(r): 



The linear-response scheme reported in Ref. J7j can be extended to all orders in 
the external perturbation. This has been carried out by Gross et al. ^1] in the general 
case of spatially inhomogeneous electron systems, and self-consistent integral equations 
for the quadratic and higher order interacting density response functions have been 
obtained by these authors. In the specific case of the uniform electron gas, which 
we are here interested in, these equations can be easily solved, to produce explicit 
interacting density response functions in terms of their noninteracting counterparts and 
the functional derivatives of the exchange-correlation potential. However, instead of 
adopting the method of solution of the above mentioned integral equations, we find 
it more instructive for our purposes, as well as self-contained, to derive an explicit 
expression for the quadratic density response function considering the uniform case 
from the very beginning. 







(2) 
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The electron density ni(r, t) + ri2(r, t) + • • • induced in an arbitrary many-electron 
system by the time-dependent external potential (p ext (r,t) coincides with the electron 
density induced in the corresponding system of noninteracting Kohn-Sham electrons 
by the time-dependent effective potential 0^(r,t) + 02^( r ?^) + " •> where 0i (r, t) is 
given by Eq. (JTJ) and 

j, e ff i 



[r,t) = J dr'v{r,r') n 2 (r',t) 
+ J dr' f dt' f xc [n ](r,t; r',t') n 2 (r',t) 
+ \jdr'J dt' J dr" J dt" g xc [n }(r,t; r' ,t';r" ,t") 



xm^OmCr",*"), (3) 

g xc [n ](r,t;r' ,t';r" ,t") being the second functional derivative of the time- dependent xc 
potential V xc [n](r,t), to be evaluated at the unperturbed static electron density n (r): 

<5 2 K c [n](r,t) 



g xc [no}(r,t;r',t';r",t"^ 



5n(r',t')5n(r",t") 

In the case of a homogeneous electron gas, there is translational invariance in all 
directions. Hence, taking Fourier transforms with respect to space and time, the exact 
momentum and frequency dependent induced electron densities, n%(q, oS) and n 2 (q, u), 
can be written as 

ni(q)=Xi(qW xt (q) = X°i(qW f (q) (5) 

and 



MQ) = J AiX2(g, qiW x \qiW x \q ~ <h) = xS(ff) 

x <tft(q) + { d% X ° 2 (q, q 1 )(f>l ff (q 1 W 1 ff (q - (6) 



where q = (q, uj), <pt (q) and 4>2 (q) are Fourier transforms of the time-dependent 
effective potentials of Eqs. (0) and 0, respectively, (p ext (q) is the Fourier transform of 
the external potential, Xi(q) an d X2{q, <?i) denote the exact linear and quadratic density- 
response functions of the interacting electron system, and Xi{q) an d X^iliQi) represent 
the corresponding density-response functions of noninteracting Kohn-Sham electrons. 
Substituting Eqs. (JTJ) and into Eqs. (J5J) and (0), one finds 

Xi{q) = r 1 {q)x\{q) (7) 

and 

X2(q, qi) = £~ l {q)xl{q, qi)^ l {qi)^ l {q - ?i) 



+ Xi(q)9xe(q, - 9i)/2, 



(8) 
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where v(q) = 47r/q 2 is the Fourier transform of the Coulomb potential, f xc (q) and 
g xc (q, qi) denote the Fourier transforms of the xc kernels of Eqs. (0) and (jlj), and e(g) is 
the test-charge-electron dielectric function |THJ UH| 

e(q) = 1 - xS(ff) Hq) + f xc {q)] • (9) 

The Fourier transforms of the linear and quadratic xc kernels in Eqs. Q and (JSJ), 
respectively, are defined as 

f xc (q) = J / TC [n ](r,*;r',Oe- <q - (r - r ' )+<t,,(t - < ' ) drA, 
9,c(q,qi) = J fe[no](r,t;r' l t / ;r ff ,Oe <Iqi - (r '- r) - Wl(t '- t)] 

x e i[( q - qiH r''-r)-( W -a, 1 )(t"-t)] dr ' ^ rf r" ^" 

Equations (J7|) and (jHJ) generalize the exaci linear density response reported in 
Ref. jEj to the realm of quadratic-response theory and the static quadratic density 
response reported in Ref. [211] to the general case of a time-dependent perturbation. 
In the so-called adiabatic LDA (ALDA), which is only rigorous in the long- wavelength 
(q — > 0) and static (uj — > 0) limits, the xc kernels f xc {q) and g xc {q, qi) are simply the 
first and second derivatives with respect to the unperturbed density of the static xc 
potential of a uniform electron gas: V^' c (no) and V" c {no). In the RPA, the xc kernels 
f xc (q) and g xc (q, qi) are set equal to zero. Finally, we find our theory in agreement 
with that of Ref. [H], while the homogeneity of the system we consider enables us to 
obtain the explicit quadratic density-response function of Eq. (JHJ) instead of presenting 
the results in the form of self-consistent integral equations. 



3. Stopping power 



There are two routes to describe the stopping power of the homogeneous electron 
gas. One is based on a perturbative expansion of the density response of the target 
(appropriate for arbitrary projectile velocities) and the other on the knowledge of the 
phase shifts of potential scattering of electrons by a statically screened impurity (only 
valid for low projectile velocities). We first consider these two alternative approaches, 
then we focus on the overlapping range of their applicability, by considering the low- 
velocity limit of the quadratic-response formulation and a second-order expansion of the 
transition-matrix elements of potential scattering. 



3.1. Quadratic density response 



To third order in the projectile charge Z-y, the average energy lost per unit length traveled 
by a recoilless probe particle moving with velocity v in a homogeneous electron gas, i.e., 
the so-called stopping power of the target is obtained as follows j^U 1221 



dE 

dx 



- 2 



zl 

7TV 



q • v 



ImXi(g) 
q 2 
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+ 



2vr 2 



qf|q- qil 2 



Im X2(g,3i) 



(10) 



which by virtue of Eqs. (0) and (jHJ) can be rigorously expressed as 



alE Z\ T r , q • v 

— = -2 Im / dq 

7Tf J 



dqi 



qfjq- qil 2 



X 



dx ir« i 1 q 2 [ q 2 ' 27T 2 

r\q) r\ gi ) r\q - gi ) [ X ° 2 (q, qi) + Xi(q)9xc(q, ?i)xS(5i)xS(g - ?i)/2] } , (ii) 
where now g = (q, q • v) and q\ = (qi, qi • v). 



3.1.1. Low-velocity limit At low frequencies we can write [H] 
and |Hj 



where 



ImXate) 9i ) = ^^q,qi + wi B qiiq + (w - u x ) B q _ qir 
0(2A; F - g) 



qi) 



.4, 



5, 



qqi 



2A r 



2iiq 

;i - fc 2 /e 



ir 1/2 



iqiiiq- qil 

x sgn (cos q )O(g/? - fci?) 
Then in the ALDA one finds 



dx irv J q 2 



A 



r 2 

q c q 



X 



Zi 



dqi 



qflqi - q| 2 



q.q^qi C q-qi ' '~ / q t q A2,(q,qi) t qi c q-qi 



^A q e q 2 XliqiXl , q _ qi /2,]}. (12) 



Here 



f q = 47r/q 2 is the Fourier transform of the Coulomb potential, hp is the Fermi 
momentum, x lq and X2,( q , qi ) denote the static (lj = 0) linear interacting and quadratic 
noninteracting density-response functions, respectively, e q is the static test-charge- 
electron dielectric function, O(x) is the Heaviside step function, q R is the radius of 
the circle circumscribing the triangle formed by the vectors q, qi, and q — qi, and 
</> q represents the angle facing q in this triangle. Evaluating some of the integrals in 
Eq. jrjj), one finds 

2k F f 



1 dE _ AZl 

V dx 37T 



dq 



<?e q 



37 + 



2Z iq 



oo rl 

dqi I dfji- 

TT JO J-l 



q qi| 2e q e qi^q-qi 



(13) 



(«q + OX2,(q, qi ) + K / ' C Xl,q 1 X?,q-q 1 /2 



2 6(g 



R 



?i|q- qil v 1 - 



■sgn(cos0 q J 
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where q and q\ now denote the magnitude of q and qi, respectively, and /i = cos0 q _ qi . 
If we put V' xc = V" c = in Eq. (fE?j) , we retrieve the RPA result of Ref . f5| . If we keep 
the actual value of V' xc but still putting V xc = 0, we reproduce the calculations reported 
in Ref. 0. 



3.2. Potential scattering 

In the low- velocity limit of a recoilless probe particle of charge Z\, the interaction 
between the Fermi gas and the probe particle can be represented as the elastic scattering 
of independent electrons by a Kohn-Sham effective static central potential V(r). Hence, 
the average energy loss per unit path length of a recoilless charged particle moving 
with velocity v (v « vp) through a uniform electron gas of density n is given by the 
following expression: 
dE 

n vk F a t r{k F ), (14) 



where 



dx 

16vr 5 r 2k 



1071" f ZK 

M*) = -£r I ^\ T fi\ 2 ( 15 ) 

is the so-called transport cross section, Tfi denoting the transition- matrix element [23J: 

T fl =< <p kf \V\i; kt > . (16) 

Here, 0k and ipk represent noninteracting and interacting electron wave functions, 
respectively, kj and k^ denote the electron momentum before and after the collision, 
k = |kj| = |k/| is the magnitude of the electron momentum, q = k/ — kj is the 
momentum transfer, and V(r) is taken to be the Kohn-Sham effective potential 

V(r) = -—+[ dv> + V xc (r), (17) 

r J | r — r | 

n(r) being the electron density induced by the presence of the static probe particle and 
V xc (r) being the xc potential at point r of the inhomogeneous electron system, which 
in the LDA is simply the xc potential of a homogeneous electron gas with electron 
density n(r). The well known expression of the transport cross-section in terms of the 
phase-shifts Si(k) of the scattering problem in the spherically symmetric potential 

4-7T °° 

°tr{k) = -7^^2(1 + 1) sin 2 ^) - 5 l+1 (k)}, (18) 
k 1=0 

greatly facilitates the numerical calculations, while the Friedel sum rule Ell 

r> oo 

Z 1 = -J2(2l + l)Si(kp) (19) 

n 1=0 

is helpful in controlling self-consistency. Echenique et al. |T21 EI] and Nagy et al. 
[T4*j evaluated the LDA Kohn-Sham effective potential by solving self-consistently the 
Kohn-Sham equation of DFT, and then computed the stopping power from Eqs. (jUj) 
and (HBJ. 
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We proceed by considering a second-order perturbative expansion of the transition 
matrix Tfi(q) entering Eq. (|15)L which will then allow us to derive Z\- and Z\- 
contributions to the stopping power from Eq. (jHJ). 

The Born series for the transition matrix element T/j in powers of the effective 
potential V(r) [which can be expanded in powers of the bare interaction —Z\jr: 
V\{f) + V 2 (r) • ■ ■] is obtained as follows 



T, 



fi 



V + VGlV + 



(20) 



where G°(r, r' 



is the noninteracting Green's function 
1 



^ik\r—r'\ 



2vr 



Up to third order in the charge Z\ of the probe particle, the square of the transition 
matrix element of Eq. (j2T)j) yields 



\T fl \ 2 = [V 1 (q)} 2 + 2V 1 ( q ) 



x 



V 2 (q)+P d qi 



^i(qi)K(q-qi) ' 
"k?/2- (1^-01)72. 



where 



and 



Z-v 



+ K + K'J ni(q) 



(21) 



(22) 



V 2 (q) = ( V(l + V^n 2 (q) 



V" r 

+ -y J dqini(qiWq- qi) 



(23) 



ni(q) and n 2 (q) being the linear and quadratic induced electron densities, and P 
in Eq. (J21)) denoting that the principal value of the integral must be taken at the 
point where the integrand is singular. It is interesting to notice that the second- 
order (Zf) contribution to |T/j| 2 has two sources. One is the first Born contribution 
to the quadratically screened effective potential V(r) and the other is the second Born 
contribution to the linearly screened effective potential, as pointed out in Ref. 0. They 
have opposite signs and it is the latter which dominates |31 IH El • 

Substituting Eqs. (0) and © into Eqs. (j22j) and and then substituting the 
expansion of Eq. (}2*Tj) into Eq. (JT5j) , one finds from Eq. (fT4*j) the following expansion for 
the stopping power: 

2k F 

j dq 
o 

. j r" ^ i:i i i:i ■) •> "n 

k\ -(k F ~ ^ < ■ (24) 
Performing the integration over the angular variables of qi, one readily reproduces 
Eq. (fTBj) . thereby proving the equivalence between the quadratic-response and potential- 
scattering schemes in the limit of low velocities of the probe particle. This generalizes 



1 dE 

v dx 



4 

3^ 
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the RPA analysis reported in Ref. |3] to the general situation where xc effects are taken 
into account. 

4. Results of numerical calculations 




0.0 1 1 1 1 1 — i 1 — = ' 

12 3 4 

r (a.u.) 



Figure 1. Stopping power of a uniform electron gas in the low- velocity limit, divided 
by the projectile velocity (friction coefficient), as a function of the electron-density 
parameter r s . The chained curves represent nonperturbative potential-scattering LDA 
calculations for protons (squares) and antiprotons (circles). The dashed (dotted) line 
represents ALDA calculations to third order in the projectile charge Z\ with (without) 
inclusion of the second derivative V" c of the xc potential. 

In Fig. ^ we plot the stopping power of a uniform electron gas of density n in the 
low- velocity limit, divided by the projectile velocity (friction coefficient), for protons 
(Zi = 1) and antiprotons {Z\ = —1) as a function of the electron-density parameter 
r s = (3/47rno) 1 ^ 3 . We evaluate both the perturbative expansion of Eq. (fT3|) and the 
non-perturbative formula§ by using the Perdew-Zunger [23] parametrization of the xc 
potential of a uniform electron gas. Our non-perturbative calculations reproduce those 
reported in Refs. ^Hj and ^1] for protons and antiprotons, respectively. Perturbative 

§ For the non-perturbative calculation, we have used the conventional scheme of iterative 

solution of Kohn-Sham equations with the potential of Eq. Ijl7(l and calculation of the stopping power 
from Eqs. i|14fl and l|18fl upon the achievement of convergence. The fulfillment of the Friedel sum rule 
of Eq. (|19fl has been monitored, the error in which has not been greater than 0.02 electrons in all the 
calculations. 
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Figure 2. Stopping power of a uniform electron gas in the low-velocity limit, divided 
by the projectile velocity (friction coefficient), as a function of the projectile charge Z\ 
and for the electron-density parameter r s — 2.07 corresponding to the average electron 
density of valence electrons in Al. The chained curve represents nonperturbative 
potential-scattering LDA calculation. The dashed (dotted) line represents our ALDA 
calculations to third order in the projectile charge Z\ with (without) inclusion of the 
second derivative V" c of the xc potential. The inset shows the same plots normalized 
to the square of the projectile's charge. 

and nonperturbative calculations are also plotted in Fig. but now for the electron- 
density parameter r s = 2.07 corresponding to valence electrons in Al and as a function 
of the projectile charge Z\. 

Also plotted in Figs. and El are perturbative calculations with no inclusion of 
the second derivative of the xc potential V xc , as reported in Ref. jH], showing that the 
inclusion of this term brings the perturbative calculations very close to the full nonlinear 
calculation in the range of high electron densities (small r s ) and small projectile charges. 
Figure [T] shows that the quadratic (perturbative) stopping power for antiprotons is 
extremely accurate for all electron densities with r s < 2. Figure 121 shows that in the case 
of Al target (r s = 2.07) and negative projectile charges the quadratic stopping power 
is accurate for the antiproton charge {Z\ = — 1) and above, but it is only accurate for 
small positive values of the projectile charge (Z% < 0.5) ||. This is due to the presence of 
the truly bound electronic states and the behavior of resonances in the case of a positive 
probe particle, which are only included in a fully nonlinear scheme. 

|| Since Z\ is the bare nucleus charge of the projectile, for non-integer Z\ the results should be 
considered as mathematical. 
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Figure 3. Stopping power of a uniform electron gas in the low- velocity limit, 
divided by the projectile velocity (friction coefficient) and as a function of the electron- 
density parameter r s , for protons (a) and antiprotons (b). Solid (dashed) curves 
represent calculations including (omitting) exchange and correlation. The chained 
curves represent the nonperturbative potential-scattering calculations. Bold (both 
solid and dashed) lines are linear (Zf) contributions. Thin (both solid and dashed) 
lines represent quadratic (Zf) calculations. 

To elucidate the role of resonances, we have performed the numerical analysis of 
the phase-shifts 5i(k) of the scattering in the self-consistent potential (|17jl. At a given 
r s with growing Z\ a resonance, which always exists in the continuum spectrum, moves 
to lower energy and grows both sharper and more intense, which results in a stronger 
variation of the phase-shifts. Since the density of states is proportional to the derivative 
of the phase-shifts the low-lying continuum states get filled preferentially resulting, 
similarly to the occupation of the bound states, in the more efficient screening of the ion 
charge and eventually in the decrease of the stopping power even before the formation of 
the bound states. On the other hand, at smaller values of Z\ the existence of weak broad 
resonances at high energies does not affect the applicability of the quadratic theory. 

In order to investigate the interplay between high-order interactions and xc effects, 
we have plotted in Fig. El the results of linear (Zf), quadratic (Zf), and fully nonlinear 
potential-scattering calculations of the stopping power of slow protons and antiprotons, 
both in the absence and in the presence of xc effects. This figure shows that: (i) The 
impact of xc effects is considerably larger within linear and quadratic response theory 
than in the more realistic nonperturbative potential-scattering approach, especially so 
in the case of protons, which indicates that there must be a large degree of cancelation 
between first, second and higher-order xc effects, (ii) At high electron densities (r s — > 0), 
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the inclusion of xc effects brings the quadratic-response calculations (thin solid lines of 
Fig. |3J) into nice agreement with their DFT based potential-scattering counterparts 
(chained solid lines). However, xc effects decrease the radius of convergence of the 
asymptotic perturbative expansion; while an unphysical negative stopping power is 
obtained within RPA for antiprotons at r s > 5.5, this unphysical behavior is obtained 
in the presence of xc effects at r s > 3.2. (iii) The performance of the perturbative 
expansion is considerably better for antiprotons than for protons. This can be attributed 
to the existence of electronic bound states and the behavior of resonances around a 
moving proton jT3|, which are out of the reach of the perturbative description, (iv) The 
performance of the quadratic response theory for the stopping power is considerably 
better than in the case of the electron density induced at the position of the projectile 
|26| I27j. which is a highly nonlinear magnitude. This is due to the fact that the stopping 
power involves an integration of the induced density over the whole space, as discussed 
in Ref. [HI- 

Finally, we note that apart from the obvious usefulness of quadratic-response 
calculations in situations where the interaction can be considered to be weak, it has 
been recently shown that perturbative calculations can be successfully used as input in 
a variational theory of charged particles interacting with a many-body system. Recent 
investigations have shown that this new variational theory brings the RPA quadratic 
stopping power for slow antiprotons into nice agreement with the corresponding 
nonperturbative potential-scattering calculations for all electron densities j28j . 

5. Summary and conclusions 

We have derived an explicit expression for the quadratic density-response function 
of a many-electron system in the framework of TDDFT, in terms of the linear 
and quadratic density-response functions of noninteracting Kohn-Sham electrons and 
functional derivatives of the time-dependent xc potential. This expression generalizes 
the rigorous linear density-response function reported in Ref. to the realm of 
quadratic-response theory, and they satisfy the self-consistent integral equations of Gross 
et al. [TT], valid for arbitrary inhomogeneous electron system. 

The exact expression for the quadratic density-response function has been used to 
obtain the stopping power of a uniform electron gas to second order in the projectile 
charge Z\, which in the low- velocity limit and within the adiabatic LDA is demonstrated 
to be equivalent to that obtained up to third order in Z\ in the framework of a 
fully nonlinear LDA potential-scattering approach. This generalizes the RPA analysis 
reported in Ref. jHj to the general situation where xc effects are taken into account. 

We have carried out LDA numerical calculations of the stopping power of a uniform 
electron gas for slow positively and negatively charged ions, as a function of both the 
electron-density parameter and the projectile charge. We find that quadratic-response 
theory yields a stopping power that is in excellent agreement with the nonperturbative 
stopping power in the range of high electron densities and small projectile charges. 
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The quadratic-response (perturbative) stopping power for antiprotons is found to be 
extremely accurate for all electron densities higher than the electron density of valence 
electrons in Al. In the case of Al, quadratic-response theory is found to yield accurate 
results for small negative projectile charges up to the antiproton charge, but a fully 
nonlinear scheme is required to account for the energy loss of slow protons. 

Although our equation (fTTj) for the stopping power is exact to the Zf order, in the 
numerical calculations we have utilized the local and adiabatic approximation for the 
linear and quadratic exchange-correlation kernels, which is consistent with the available 
fully nonlinear calculations within the potential scattering method. To study the role 
of the non-locality (wave-vector dependence of the exchange-correlation potential) and 
non-adiabaticity (its frequency dependence) in the nonlinear theory of stopping-power 
is, however, a challenging task, and this work is now in progress |29j . 
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